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I. INTRODUCTION 

Several authors have noticed that the free space 
Maxwell equations in an arbitrary co-ordinate system, 
on a Riemannian space-time background, take the same 
form as they do in a certain class of transparent, inho- 
mogeneous, anisotropic media (see references in [JJ and 
also [2 §90]). Recently this relationship has been put to 
work in reverse. Transformation optics [3H7] has devel- 
oped to use this formal analogy to design optical devices, 
and even investigate analogues of astrophysical objects 
in the laboratory [8j |9] . 

Light propagation on a Riemannian space-time back- 
ground can be intuitively understood in terms of rays 
following geodesies, and polarization undergoing parallel 
transport along each ray [3J [10] . Transformation optics 
uses this simple picture to design a space that acts on 
the optical field in a desired way, and, using the peculiar 
property of Maxwell's equations mentioned in the pre- 
ceding paragraph, thereby determines the necessary ma- 
terial properties from geometrical quantities. Notably, 
this recipe has been applied to derive the material prop- 
erties necessary for devices that conceal objects from the 
electromagnetic field OUT], and focus light intensity into 
a region that is smaller than the diffraction limit leads us 
to expect [TU [T3] . In this paper, the idea is to provide 
more variables for this design strategy to explore: in par- 
ticular, to incorporate chiral media into transformation 
optics. 

As stated above, the existing theory of transforma- 
tion optics works within Riemannian geometry, where 
the 'design parameter' is the space time metric, g^ u . 
Yet, the term transformation optics came from the initial 
use of co-ordinate transformations to arrive at material 
parameters — i.e. Euclidean geometry. Therefore, in the 
original sense, transformation optics works through the 
specification of three functions of position. Clearly the 
full Riemannian geometry has greater freedom, with a 
symmetric space-time metric containing ten independent 
functions of position, translating into nine independent 
material parameters [35] , 
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However, even if we recognise that a geometry must 
affect electric and magnetic fields in the same way, 
then this may not be the full story. Symmetric, 
impedance matched permittivity and permeability ten- 
sors (i.e. e/e = /V^o) represent six independent quan- 
tities, and the possibility of magneto-electric coupling, at 
least another six components. This counting argument 
leaves three real magneto-electric coupling parameters 
that cannot be represented within Riemannian geometry, 
but may well have a geometrical interpretation |44j . It 
is therefore worth investigating non-Riemannian exten- 
sions to transformation optics. In such a geometry there 
is at least one additional field — the space-time torsion, 
T^ va — T^ vrJ — r^o-y — that may be freely specified inde- 
pendent of the metric. Torsion has already been explored 
in other analogue systems: e.g. the theory of sound waves 
propagating through superfluids, where non-Riemannian 
geometry has been used to describe the interaction with 
vorticity [TBI [17] . 

In the following it will be shown that if we couple 
a non-Riemannian geometrical background to the free 
space Maxwell equations in a certain way, then Maxwell's 
equations can be interpreted as if in an inhomogeneous, 
isotropic, chiral medium described both by a Tellegen pa- 
rameter, x an d a chirality parameter, k. In the limit of 
geometrical optics, this coupling is shown to reproduce 
the usual geodesic and parallel transport equations, but 
in the presence of geometrical torsion: i.e. optical activ- 
ity is shown to have a geometrical interpretation in terms 
of the torsion tensor. 



II. THE RELATIONSHIP BETWEEN EM 
FIELDS IN CONTINUOUS MEDIA AND 
RIEMANNIAN GEOMETRY 

We begin by reviewing the existing theory of transfor- 
mation optics [ni m m] • 

Consider a space-time that is not necessarily flat — i.e. 
where the curvature tensor, i? M 1 y (TT , may not vanish — 
and where the co-ordinates, x^, are arbitrary. In this 
instance, we write down the free space Maxwell's equa- 
tions using their usual four dimensional form [2\, along 
with the convention from general relativity |19j . that or- 
dinary partial derivatives be replaced by covariant ones, 
dfj, — > V M , and that the permutation symbol be scaled by 
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the volume element, e^ vaT -> e ^ar - 

e^° T V v F aT = 
V^F"" = 0. 



-1/2 P.VOT 



(1) 

(2) 



The covariant derivative, V„, in dll [2[) differs from an or 



dinary partial derivative by a quantity, T a vfJi , known as 
the connection symbol [5U] . Transformation optics works 
because, in the Riemannian case, the connection symbol 
in ([l}j2| plays the same algebraic role as the difference be- 
tween the field equations in vacuum and in a polarizable 
medium. We can see this through explicitly identifying 



rV in Q and 
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r V Q = o. (4) 



The assumptions of Riemannian geometry lead to a par- 
ticular form for that is known as a Christoffel sym- 
bol [20], r^i^cr = an d depends only on the form of 
the metric tensor, g^, 
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Due to the symmetry of the Christoffel symbol in its 
lower two indices, ^ is indifferent to the distinction be- 
tween and d^. Therefore the definition of the field 
tensor in terms of the vector potential is the same as if 
the co-ordinates were those of a Galilean system [45] . 
Furthermore, the antisymmetry of the field tensor means 
that the final term to the left of the equals sign in Q 
is also zero. We can now see that the only term that 
distinguishes the Maxwell equations ^ and Q from 
a Galilean system is proportional to the trace of the 
Christoffel symbol g0], 

where g = Aet{g^ v ). The free space Maxwell equations on 
a Riemannian background, and (|2|, can thus also be 
written in a very similar form to the free space Maxwell 
equations in Galilean space-time. All that is changed is 
the appearance of a factor of y/—g~, and the relationship 
between F^ u and F„„ 06], 



e^ T d„F aT = 



(6) 



-gF^) = 0. (7) 

There are two equivalent ways to understand ^ and 
([7]) — either as we have done so far, in terms of an empty, 
possibly non-flat, space-time background; or equiva- 
lently in terms of a Galilean system, containing a dielec- 
tric medium. For, if we consider the four dimensional 
Maxwell equations in the presence of a material medium, 
within a Galilean co-ordinate system, and without any 
external sources, then Maxwell's equations can be writ- 
ten as |2"51. 
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In this case there are two separate, but not independent, 
'field tensors'; the tensor, F aT — a bivector containing the 
physical fields, F aT = (E/c, B); and the tensor, — a 
bivector containing the material fields, G^ = (— cD, H), 
the vanishing divergence of which indicates that the ma- 
terial has no net charge, and that no net current passes 
through any cross section. The constitutive relationship 
between F aT and G^ T serves to indicate how the physical 
fields are influenced by the medium. 

As the form of ([6]-[7]), is identical to that of ([8jj9]), there 
is a correspondence between each and every co-ordinate 
system on any Riemannian background and an equivalent 
material, described within a Galilean co-ordinate system, 
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where the factor of /Iq has been introduced so that 
the units agree with the physical interpretation of G^" . 
The co-ordinates, x M , on a general background, are re- 
interpreted as a Galilean system, g^ v f-> 77^, with the 
contravariant field tensor appearing as a material field. 
The constitutive relationship between the physical fields, 
F^, and the material fields, G^ LV is obtained from the 
relationship between covariant and contravariant indices 
on the Riemannian background, 



F^iv 9[1<79vtF 
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G aT , 



(11) 



To make the interpretation more transparent, we intro- 
duce the usual three dimensional quantities; Ei = Fq;; 



B, 



» k F jk ; cD z = G M ; and H t 



adopt dyadic notation, so that (11) becomes 



D = e E 
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B = (i-H- lg xE 



where, e, 



G^, and 

(12) 
(13) 

and we 
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Hi = (9oi9oj/9oo ~ 9ij)\ 
The conclusion is there- 



have introduced the symbols; 

lij 1 = ~9 %J 'i and g 3 = cg oj /g 00 - 
fore that each co-ordinate system, on every Riemannian 
background, can be considered to appear to the electro- 
magnetic field as a continuous medium with equal rela- 
tive permeability and permittivity tensors (an impedance 
matched medium) — eij/eo — Mij/Vo — described within a 
Galilean system of co-ordinates. 

Perhaps even more importantly, the reverse also holds; 
for a fixed frequency of the electromagnetic field, all 
transparent, impedance matched media can be under- 
stood in terms of the vacuum Maxwell equations within 
a Riemannian geometry. However, this space-time co- 
ordinate system is not uniquely defined by the medium, 
due to the invariance of the Maxwell equations under 
conformal transformations. 



3 



III. ISOTROPIC CHIRAL MEDIA 



Section IV will describe a class of non-Riemannian ge- 
ometries that are equivalent to inhomogeneous, isotropic, 
chiral media. Unfortunately, there seems to be no agree- 
ment on the form of the constitutive relations that should 
be used to describe such media. Although, in the fre- 
quency domain, these various constitutive relations can 
be shown to be physically equivalent [241 125) , the mean- 
ing of the individual material parameters is different for 
each constitutive relation. This will turn out to be im- 
portant when we come to interpret chiral parameters in 
terms of geometrical quantities. Therefore we spend this 
section distinguishing the constitutive relations before, in 
the next, introducing the geometry. 

Landau and Lifshitz determine the relationship be- 
tween the material fields and the physical fields in op- 
tically active media in terms of anti-symmetric complex 
components in the permittivity tensor: = e* t [2"3"] . 
This is equivalent to what is known as the Drude-Born 
constitutive relation [2H [5S], where D is coupled not 
only to E, but also to V x E. For a geometric inter- 
pretation, we must have constitutive relationships that 
are local and have symmetric coupling terms in D and 
B (or H) [37], and so exclude this possibility. Instead, 
we describe isotropic chirality via the electric-magnetic 
coupling terms that were discussed in the introduction. 
The Tellegen constitutive relations are often used [24 , 
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B = /x T ■ H + - (xt + in-r) E. 
c 



(14) 
(15) 



where x represents the Tellegen parameter [57] , and k the 
chiral parameter. The presence of the i indicates that, 
so long as K 7^ 0, (14 1 and (15 1 only have a meaning in 



the frequency domain. While the existence of media with 
non-zero k is unquestionable, there is a history of debate 
regarding the reality of media with non-zero x [55H3"T] . 

We can observe that an isotropic chiral medium is a 
material where there is a linear coupling between one 
component of the electric (magnetic) polarization, say, 
P x (M x ), and the same component of the magnetic (elec- 
tric) field. Therefore, for our purpose it is useful to see 
our initial suspicions confirmed, and to notice is that it is 
not possible to have a totally spatially isotropic medium, 
and have such a direct coupling described geometrically 
by (fT2| and (13), for in this case g should vanish. Hence, 



optical activity in isotropic media cannot be understood 
in terms of the theory of transformation optics presented 
in section |TT] 

Part of the ambiguity in the form of the constitutive 
relations for these media comes in deciding whether to 
use B or H to describe the coupling of the magnetic 
field to the electric polarization. One common alternative 



form to (14 151 is the Boys-Post relation 
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In the frequency domain (where these relationships are 
defined) (14 15) and (16 17) are equivalent. However, it 



is important that in each case the meaning of the per- 
mittivity and permeability is different, as well as the in- 
terpretation of the chirality. Indeed, if we cast ( 16 1 and 



(17) into the form of (14) and (15 1, then we obtain 
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E + - (xp + inp) /x P ■ H 
c 



/Lt p • H + - (xp - iK P ) /x F 
c 
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so that we can observe a correspondence in the coupling 
parameters, 



e T e P + 4f (xp 

/x T <-» /x p 
Xt ± in T <-» /x p (xp ± i/Sp) 
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Therefore, in general, when the magnetic susceptibility is 
anisotropic, whether the chirality is isotropic is relative to 
the interpretation. Furthermore, 'impedance matching' 
in the Post constitutive relations does not translate into 
impedance matching in the Tellegen interpretation. As 
noted in |24j . true impedance matching must be done 
relative to the Tellegen parameters. However, in the limit 
of small chirality, the difference is negligible. 



IV. NON-RIEMANNIAN GEOMETRY AND 
ISOTROPIC CHIRAL MEDIA 

From the summary of the theory of transformation op- 
tics given in section[H] it is evident that there are two rea- 
sons why a Riemannian background has the same effect 
on the Maxwell equations as a material medium; (i) the 
definition of the field tensor in terms of the vector po- 
tential is unchanged by the background, so that the co- 
ordinates may be re-interpreted as if they were Galilean, 
and; (ii) the expression for the trace of the Christoffel 
symbols is such that Q can be written as a total di- 
vergence. Consequently, when the co-ordinate system is 
re-interpreted as Galilean, the effect of the geometry is 
that of a charge neutral medium, with no net current 
passing through any cross section. 

Is it possible to fulfil both these conditions with a more 
general geometrical background? Let us explore more 
general forms of the connection, 



{la} + C^vo 



(19) 



It is instructive to write the non-Riemannian part of the 
connection, C^w, in terms of component tensors (e.g. 
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FIG. 1: A completely antisymmetric torsion field, 



-Tvtio, does not affect geodesies, and rotates vectors 



in the plane normal to the direction of parallel transport. In 
the figure, parallel transport is being performed along the blue 
lines. The transparent vectors show the case when torsion is 
not included, while the solid vectors show the case when we 
have a completely antisymmetric torsion field. 



see [IH1 Section 3.3]), 



where K^ vrr = 



[T av <j + T uaa + T aav ] is the contor- 
sion tensor (recalling that the torsion, T^ VG is the anti- 
symmetric part of the connection), and H^ ua represents 
the non-metricity [35]. When H^ ua = 0, the connection 



satisfies the condition, V ^g aT = 
follows we assume that H^ va = 



0, everywhere. In what 
, so that, C» = K»„„. 



Almost immediately we hit problems if we add „ CT 
into the connection. For if we follow the usual 'par- 
tial derivative goes to covariant derivative' rule, then the 
anti-symmetry in the lower indices interferes with the 
definition of the field tensor, ([6]), 



iivor rg p _pp pi 



0. 



In itself this might not be a problem, did it not break 
the gauge invariance of the theory. For if we perform a 
gauge transformation, A'^ = + V M y>, the field tensor 
ends up depending upon ip, 
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Therefore, we reach the conclusion that, with non-zero 
torsion, the 'partial derivative goes to covariant deriva- 
tive' rule does not produce a gauge invariant theory. For 
more extensive coverage of this issue, see [HI Section 
11.3] and [33] , Hence we choose to keep the definition 
of the field tensor the same as in the Riemannian case, 
Fpu = dpA v — d v Ap. 

The second Maxwell equation does not relate to the 
definition of the field tensor, and we may suppose that the 
background geometry modifies this equation with terms 
including the contorsion tensor. Applying (19) to Q, we 



obtain. 



^=d^ (V^gF^) + \t v ^f to = o. (20) 
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In order that this appear as a material medium — c.f. 
([9]) — there is a minimal choice for the form of the tor- 
sion, 



(21) 



where x is an arbitrary single- valued function of x M . 
The Maxwell equation associated with sources, ( 20 ), thus 
takes the following form, 



-gF^ + \xw v ° T »F TC 



(22) 



It is interesting that ( 22 ) has the same form of coupling 



to the torsion field as was obtained in [33 through micro- 
scopic considerations of the interaction between a space- 
time torsion field and the vacuum polarization associated 
with the quantized electromagnetic field (see figure [I] for 
the geometrical interpretation of (21)). The divergence- 
less material field in (22) is, 



Mo 



+ X - X e^F ar (23) 



Due to the appearance of the dual electromagnetic field 
tensor on the right hand side of ( 23 ) , the presence of the 



torsion (21 ) in the connection is equivalent to some cou- 



pling between like components of the polarization (mag- 
netization) and the magnetic (electric) field. To be ex- 
plicit, we write (23) in terms of the usual fields, (E, B, 
D, H). 

Take the simplest case first, and consider an isotropic 
material, g^ v — diag(l, — n 2 , — n 2 , — n 2 ). For this case 
(23) becomes, 
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(24) 



The addition of the torsion, (21), into the isotropic 



geometry may be understood in material terms as an 
isotropic chiral medium with a non-zero Tellegen pa- 
rameter, that has equal relative permittivity and perme- 



ability as regards the Boys-Post prescription, (16-17), 

/i p /jUo = ep/eo- 

The next simplest case is where the medium is 
anisotropic, but where the magneto-electric coupling de- 
fined by the g vector vanishes. Here, (23) gives, 



D=e -E+-B 
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(25) 
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where, as before, ti 



ill 



and jXij 
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As the torsion, (21 ), is actually a pseudo-tensor, the \ pa- 
rameter changes sign under time reversal or spatial par- 
ity inversion: this is consistent with a Tellegen medium, 
which is not time reversible. The constitutive relation is 
somewhat more complicated in the case of a space-time 
background with a non-vanishing g vector as there is an 
interplay between the magneto-electric coupling due to 
the metric, coming in the form of the g vector, and that 
coming from the torsion. 

We conclude that the inhomogeneity of the Tellegen 
parameter is equivalent to space time torsion in the co- 
variant derivative of Q. Notice that when the Tellegen 
parameter is homogeneous, it disappears from the geom- 
etry, which again becomes Riemannian: this is consistent 
with the known invariance of the Maxwell equations un- 
der transformations of the fields, where, D — > D + rjH 
and H — > H — 77E, with uniform r\ (e.g. see [2"51 151] ). 

It is evident that the more physically important 
magneto-electric parameter, k, does not arise from the 
above modification to the space-time connection. The 
definition of this quantity, (16-17), anticipates that it 
can only arise in the frequency domain, which we now 
consider. 



A. The frequency domain 



From the form of ( 28 ) and ( 29 1 , the torsion given by ( 27 ) 



defines the following divergence-less quantity, 




In terms of a Galilean system containing a material 
medium, these are equivalent to the vector relationships, 
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e E + i-B 
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t 1 



■B + i-E, 
c 



which define a material medium with a chiral parameter, 
k, interpreted in the sense of a Boys-Post constitutive 
relationship. Combining this result with that of the pre- 
vious section, a medium with both a chiral parameter 
and a Tellegen parameter may be defined via a torsion 
pseudo-tensor with the components, 



rpijk 


0,, 

c 
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= w ijk dk (x ~ 
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= Moe^d* ( X + 


in) 



(30) 



In the frequency domain, we consider a purely spatial 
geometry — goo = 1, got = gia = — and replace the time 
derivative with — ito. With this assumption, (201 is, 



(26) 



where we have assumed, as in the previous section, that 
the additional torsion does not alter the trace of the con- 
nection symbol. One such set of components are as fol- 
lows, 



As we are in the frequency domain, we have assumed 
that the material parameters are independent of time. 
Notice that when the material is uniform, the mixed time 
and space components of the T"^ 1 vanish, and we are 
left with only the spatial components of the torsion that 
are proportional to k. Indeed, the spatial torsion is by 
far and a way the dominant part of the object, as it is 
also weighted by the factor, u/c. This is particularly 
important in the limit of geometrical optics, which we 
must now investigate. 
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Note that, despite appearances, (27) transforms as a 



pseudo-tensor under purely spatial co-ordinate transfor- 



mations. We can decompose (26) into two equations, one 
for v = 



0, 



/'■o 







(28) 



and one for v = j, 
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V. PARALLEL TRANSPORT AND 
GEOMETRICAL OPTICS 

Transformation optics is a geometrical theory that goes 
beyond ordinary geometrical optics: as illustrated above, 
it is concerned with an exact mapping of Maxwell's equa- 
tions from a geometrical background onto an equiva- 
lent material medium. However, geometrical optics en- 
codes a great deal of the intuitive content of the theory: 
rays follow geodesies, and polarization is parallel trans- 
ported ^] . Therefore it is a minimal requirement that 
the theory of section \FV\ bear these intuitions out: geo- 
metrical optics should behave as a theory of rays on a 
background with torsion. Here we show that this is in- 
deed the case. The approach of this section is like that of 
the last, the existing theory is briefly reviewed so that we 
can then bring out the features of the non-Riemannian 
modifications. 
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A. Riemannian media 

The starting point of geometrical optics is the wave 
equation, which in a non-chiral, impedance matched 
medium, can be obtained from the Riemannian form of 
Maxwell's equations, ([6[[7|. As noted previously, it is im- 
material whether covariant or partial derivatives appear 
in the definition of the field tensor in a Riemannian ge- 
ometry. Therefore the four-dimensional curl of ^ can 
be written as, 



(31) 



As has been assumed throughout, the covariant deriva- 
tive of the Levi-Civita symbol is zero, and the usual for- 



mula, 6*22 



c ai3-yfJ, 



-where 5"27 is a 3 x 3 deter- 



minant of Kronecker deltas [30] — can be used to obtain 
a sum of second derivatives of the field tensor, 

V 7 V"F' 37 + V 7 V 7 F Q ^ + V 7 V /3 F 7 " = 0. (32) 

Applying ^ then gives the wave equation, 

V 7 V 7 F Q(3 + [V 7 , VJ F Pl + [V 7 , Vp] F ia = 0, (33) 

where the commutators of the derivatives are propor- 
tional to the contraction of the Ricmann curvature ten- 
sor, W vr , T against each index of the field tensor in turn: 
e.g. for a four-vector, V 1 , [V a ,V^]Vo- = R T apaYr- 
These curvature terms will not enter the approximation 
of geometrical optics, because they do not diverge as the 
wave-length goes to zero. To see this we write down 
the field tensor in the form, F^ = / M „e 27rl,s / A , where 
A = 2-kc/lj, and take the limit of rapidly varying phase: 
A — > 0, which is equivalent to assuming a length scale 
for the variation of the material properties that is much 
larger than the wavelength of the optical field. 

Expanding (33) in powers of A -1 , and requiring that 



the coefficient of each power vanish separately (neglect- 
ing the zeroth order as A — > 0), yields the equations of 
geometrical optics, 



A" 



(d a s) (d a s) = o 



(34) 



A" 1 : (d 7 S)Vy^ = -/^V 7 (d 7 S) (35) 



which, as expected, do not include the curvature tensor. 



The first of these relations, (34) is equivalent to the state- 



ment that rays follow geodesies, for if we operate on the 
left with the covariant derivative, V^, and remember that 
in a Riemannian space-time, the second order derivatives 
of a scalar commute (the torsion is zero): [V Q , V/3] <p — 0, 
then, 



(d a S) S7p (d a S) 



D 2 xp 
ds 2 



= 0, 



where the gradient of the phase, d a S is taken to equal 
the tangent vector to a curve, d a S = g a pdx@ /ds, and 
D = dx a W a . Equation (36) is precisely the rule for the 



geodesic motion of a material particle in a Riemannian 
space-time. 

The second defining equation of geometrical optics, 



(35 ), illustrates how the polarization changes along a ray. 



We write the field tensor amplitude as a bivector that, 
in the frame we are considering, contains two unit three- 



vectors, multiplied by an amplitude, /„,, 



(u, 



UiU" 



T, where, 
0. 



1; ViV 1 = 1; Uik 1 = V{k % = U{V 
The tensor, indicates the direction of the polariza- 
tion, and p5|) becomes, 



dx J \ 



2 / Du^ 
I ds 



= 0. 



(37) 



Both terms to the left of the equality in ( |37| ) must vanish 
separately if the unit polarization vectors defined within 
are to remain unit vectors at all points along a ray 
(i.e. otherwise they would grow or diminish exponen- 
tially with the relative change in the energy density). 
Consequently, 



ds ) 
Du,,,, 



ds 



(38) 
(39) 



Explicitly, ( 38 ) is equivalent to the continuity of the en- 



ergy momentum along the ray, and ( 39 1 to the statement 



that the direction of the polarization is parallel trans- 
ported along the ray. These are the essential results of 
geometrical optics in a Riemannian medium. 



B. non— Riemannian media 

If we allow for the possibility of a space time with tor- 
sion as discussed above, then it is not immediately obvi- 
ous whether the equations of geometrical optics will still 
carry the simple structure where rays follow geodesies, 
and polarization is parallel transported: although the 
Maxwell equations remain identical in form to ([6]) and 
([7]) throughout section IV this was achieved through im- 
posing the definition of the field tensor, and we should 
check that this imposition has not spoilt the geometrical 
interpretation. We shall show that, due to the particular 
form of the torsion, (30), the geometrical optics limit is 



unaffected; rays still follow geodesies, and the polariza- 
tion is still parallel transported. 

The analogous situation to that represented by equa- 
tion (31) contains an additional contribution that arises 



because of the use of partial, rather than covariant deriva- 
tives to define the field tensor, 



T a va F aT ) = 



Following the same procedure that led to (32), and ex- 



(36) panding the contraction of the Levi-Civita symbols 



V 7 V7 F «/3 + [y 7i ya] pPj + [y^ y/3] pi* 

V 7 (r p a ~<F pl3 + T/ a F pl + T p 7/3 F pQ ) = 0, (40) 
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it is clear that there are terms in addition to the wave 
operator, V 7 V 7 F Q ^, that involve not just the contrac- 
tion of geometric quantities against the field tensor, but 
also against derivatives of the field tensor. Such terms 
will remain in the equations of geometrical optics, and 
we must investigate them further. 

The commutator of the derivatives of the field tensor, 
[V 7 , V a ] F^ , now also contains terms involving the tor- 
sion tensor, 

[V Q , V/,] FK» = W aaf) F° v + R\ ap F^ + TP a0 V p F^. 

(41) 

Using (41| and applying Q, we can group the terms in 
addition to the wave operator in (40). The zeroth order 



the meaning of these equations is now slightly different. 
For if we take (34) and attempt to derive the geodesic 



equation as before then wc find an additional term, 



(d a S)V p (d a S) 



(d a S) V a (d p S) + (d a S) [V/3, V a ] s 
(d a S) V a (dpS) + T al3 ^S)(VS) 



For the same reason that this limit works in the first 
place, namely the vanishing of the terms in (42), this 
additional contribution vanishes and we have, 



ds 2 



ds 2 



dx a dx a 
ds ds 







(43) 



terms involving no derivatives of the fields are found to 
be, 

th : g aa (BP^F" + Wp^P) 

- g Pa {B?„ F" + W pia F a ?) 

+ F p/3 V 7 T p Q7 - F pQ V 7 T p 7 ' 3 + F^V 7 T/ a 

meanwhile, the terms that are first order in the deriva- 
tives of the field tensor are, 



This is the equation for an auto-parallel rather than a 
geodesic Section 10], as it contains the full connection 
and not only the Christoffel symbol: this equation for- 
mally determines the straightest line between two points 
and not the shortest. However, again due to the anti- 



symmetry of (30) in all indices, the contribution of the 
contorsion is zero (i.e. chiral media are equivalent to a 
geometry with a K^va that is antisymmetric in the lower 
two indices), and rays follow geodesies: r /3 ' aa — > {£ a \. 



The derivation of the equivalent of ( 37 1 is unaltered in 



this situation, and so (38) and (39) remain in the same 



I s 



T w/5 ) V 7 F p a . 

(42) 



form. Firstly the propagation of energy-momentum, 



T 



> dx 1 
~ds~ 



> dx 1 
ds 



dx a 
'~ds~ 



0. 



For a general form of T"^ 7 , (42) is non-zero, and a sim- 
ple geometric description will not apply. However it is 
immediately clear that the description of Tellegen media 
given by (21 1 makes these first order terms vanish. 



The contorsion that gave rise to the chirality did not 
alter the trace of the connection. Therefore energy- 
momentum propagates relative to geodesies, as in the 
Riemannian case: this is consistent with the equivalence 
Intreating (|42] in the frequency domain, we assume, as of ^ tQ geodesic motion . Meanwhile, the propagation 



m [6&\ . that the optical activity of the chiral parameter, 
k is such that the polarization is only slightly changed 
over each optical cycle: i.e. that ujk/c does not diverge 
as A — > 0. This means that the quantity k is of order 



of the polarization along the ray is affected by the pres- 
ence of the torsion, 



Du 



du 



pu 



A, and the non-zero part of (30) in the limit equals only 
the spatial part of the torsion, (not including the 

Tellegen parameter, the contribution of which we have 
shown to equal zero), so that (42) also vanishes in this 
case. This proves that the limit of geometrical optics 
involves only the wave operator and zeroth order terms, 
just as in the Riemannian case. 

We have therefore established that, as in the usual sit- 



ds 



ds 



pa U'ov 



dx a 
ds 



r r 



VOL "^(7 



dx a 
ds 



0. (44) 



uation presented in section VA only the wave operator 
matters in the limit of geometrical optics, 

V 7 V 7 F^ = 0. 

For slowly varying torsion and small curvature and tor- 
sion in comparison to 1/A this is the equation obeyed 
by the exact solution to Maxwell's equations, to a good 
approximation. The geometrical understanding of the 
theory becomes more complicated when the torsion is 
rapidly varying, just as it does when the curvature is 
large in the usual theory of transformation optics. So 
our formalism passes the first test. 

Inserting the ansatz for the field tensor as in the pre- 
vious section, we have, again (34) and (35). However, 



So the formalism appears to be consistent with the idea 
that weak chirality should act only to rotate polarization 
during propagation (c.f. figure [T]). 

In summary: for a chiral medium that rotates the po- 
larization by a finite amount over a typical length scale, 
and where the change in the chiral parameter, k, is not 
significant over a wave-length, the geometrical optics of 
chiral media requires that we add torsion into the con- 
nection. The form of the torsion, (30), is such that the 
propagation of a ray is unaffected — geodesies are equiv- 
alent to auto-parallels — while the parallel transport of 
the polarization is modified. This result also holds for 
Tellegen media, where \ does not vary too rapidly. 



VI. APPLICATIONS 
A. A homogeneous, isotropic, chiral medium 

The simplest test of this theory is to apply it to the 
simplest kind of chiral medium; one that is homogenous. 
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Note that throughout this section and the next we im- 
plicitly work in the frequency domain, despite using four 
dimensional notation. For an isotropic, homogeneous 
medium, the theory of sections |IT| and [TV| prescribe that 
the metric and torsion should be given by, 



9^ 



1 













— n 2 













-n 2 













— n 



and, 



T 0ij = T tOj = q. T ijk = _ 2m 



o — e 
c 



ijk 



Due to the assumed uniformity of n, the Christoffel sym- 
bols vanish, and the connection is equal to the contorsion 
tensor, 



Therefore (43) becomes 



d 2 x? 
ds 2 



= 



(45) 

Rays follow straight lines. Yet, the polarization is 
changed along each ray, as is clear from (44 1. For in- 
stance if we take the unit vector for the electric field, 
uoi = Ui, then it changes according to, 



du,j 
ds 



(46) 



where a = hqluk/c and kj — dxj/ds. Suppose that a 
ray travels along the x-axis, k = (n, 0,0). In this case, 
u = (0, u y , u z ), and (46) yields two coupled equations, 



du y 

ds 
du z 

ds 



(47) 



From the definition of dx l /ds, dx 1 dxi — ds 2 , the line 
element on the ray is, ds = ndx. Therefore the covariant 
unit vector, u has the following form, 

u = n (0, sin (nax), cos (nax)) . (48) 

An almost identical calculation for the unit vector of the 



magnetic field, v shows that v also satisfies (47). Apply- 
ing the definition of the field tensor, 



n (0, — cos (nax), sin (nax)) 



(49) 



From ( 45 1 , ( 48 1 , and H9h , the field in the material is 



proportional to, 

E = (0, sin (nax), cos (nax)) e^ {nx - ct) (50) 
(0, - cos (nax), sin (nax)) e ^ {nx - ct) (51) 



B 



n 
c 



In the weakly chiral limit in which we are working, (50) 
and ( p3l] ) are the solutions to the Maxwell equations. To 
see this we consider the wave equation that arises from 
the usual Maxwell equations (/i/ u o = e/ e o) with the con- 
stitutive relations (16) and (17), 



V 2 E 



V x E 



-E = 0. 



(52) 



Substituting in an electric field of the form u(x) e 1 
into (52) gives, 



du iiQnujK 

— + xxu 

ox c 



ic 
2nuj 



d 2 u 

dx 2 



2uQnnui 



V x u 



(53) 



The right hand side of ( 53 ) is proportional to A times a 
quantity of order unity. Therefore, in the approximation 
of (l50T[5ll, 



du 

dx 



fiQnuK 



0. 



which is identical to ( 47 ) . This proves that in the case of 
an isotropic, homogeneous, chiral medium, with uqk <C 
1, non-Riemannian geometrical optics is equivalent to 
the solution of Maxwell's equations. 



B. Maxwell's fish eye lens 

In the formalism of transformation optics presented in 
section [TlJ geometry is implemented for the purpose of 
directing rays, and polarization is a bystander, and must 
respond in a way that is determined by the geodesies. 
However, one may wish to maintain a given polariza- 
tion throughout a device, or change it in some prespec- 
ified manner. Here we show, using the simplest exam- 
ple of a curved geometry for light — the Maxwell fish 
eye [TOl H3 [35] — that torsion can be used to control 
the polarization of light without affecting the geodesies. 
We should note that chiral media have previously been 
considered for 'correcting' polarization in a variant of the 
planar fish eye [33] . 

The fish eye is a continuous medium in which the be- 
haviour of light maps onto the free motion on the surface 
of a sphere (this may be a 2-sphere or a 3-sphere, de- 
pending upon whether the medium is planar, or truly 
three dimensional). It has recently recieved much atten- 
tion, in part due to its ability to periodically perfectly re- 
construct an initial optical pulse as it propagates through 
the medium, as well as sub-diffraction resolution (T3l I39T 

If the behaviour of light in terms of lab co-ordinates 
(x,y,z) corresponds to motion on the surface of a 3- 
sphere (a hypersphere) , the optical line element can be 
written in a form corresponding to an isotropic medium, 



dV 



n(r) 2 (dx 2 
2 [dQ 2 + 



sm 



dy 2 + dz 2 ) 
2 (9) (d® 2 



! (d>)d x 2 )] , (54) 
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If we introduce spherical polar co-ordinates in the lab 



(a) 





(b) 

FIG. 2: In Maxwell's fish eye lens, the ray trajectories are the 
great circles of a 3-sphere, as illustrated for the special case 
of r — a in (a) and (b). (a) shows the usual situation, where 
the refractive index determines both the ray trajectory (red 
line) and the behaviour of the polarization (blue lines). Mean- 
while (b) shows that transfo rma tion optics with the inclusion 
of the geometrical torsion, (58 1 (i.e. isotropic chirality), al- 



lows for the manipulation of polarization without affecting 
the geodesies. 



where r = \/ x 2 + y 2 + z 2 , and, a, (0,$,%), are the ra- 
dius and angular co-ordinates of the equivalent 3-sphere. 

A spherical geometry has no boundary, whereas phys- 
ical space has a 'boundary' at infinity. Therefore, as 
r — > oo in physical space, n(r) should be such that the ra- 
dius of any circle surrounding the origin, as experienced 
by a light ray, 2mi(r)r, goes to zero. Also, the length of 
any optical path, J Q n(r)dr, should equal ira, where a 
is the radius of the sphere. A refractive index fulfilling 
both these requirements is, 



n(r) 



2a— arctan (r la) 
dr 



1 + {r/af 



(55) 



system — dx + dy 2 + dz 2 



dr 2 



\d0 2 



■(O)d^) 



in (54) — then the observation that the angles, 8, (j), must, 
by symmetry, equal the corresponding angles on the 3- 
sphere, 9, (j> — 3>,x leads, via (55), to the identification: 
r = atan(0/2). Performing this transformation of the 



radial co-ordinate in the first line of (54) gives the line 



element on the second line of (54), justifying (55) 



The metric tensor associated with Maxwell's fish eye 



can be immediately written down from ( 54 ) , 



9nv 



(I 

-a 2 



\0 



-a 






2 sin 2 (8) 




; sin 2 (e)sin 2 ($)/ 



This corresponds to the non-zero Christoffel symbols; 
{22} = -sin(e)cos(9); {^g} = - sin(6) cos(6) sin 2 ($); 
{ 2 2 } - cot(6); { 2 3 } = -sin(a>)cos($); {^} = cot(6); 
and {23} = cot($). These symbols determine the tra- 
jectories of the light rays, as well as the change in the 
polarization during propagation. In particular, light rays 
with motion along the <!>, x axes are confined to a 2- 
sphere (of constant 0) when <d = ir/2 (r = a in the lab 
system). As a visual example, let's examine the motion 
of light rays and transport of polarization in this specific 
case (see figure [2|. 



The geodesic equations (36) on the 2-sphere defined 
by 9 = 7r/2 are, 



d 2 Q 
ds 2 
d 2 $ 
Is 2 

<Px 

ds 2 





sin ($) cos ($ 
-2 cot ($) 



dxdx 
ds ds 
d<& dx 
ds ds ' 



(56) 



so the rays remain confined on the 2-sphere. The change 
of the polarization is similarly obtained, through apply- 
ing pi 



du 



01 



ds 
du 02 

ds 
du 03 

ds 



= 



dx 



= COt ($)m 3 , 

ds 

d$ dy 

= cot ($)u 03 - J sin ($) cos ($)m 02 ^- 

ds ds 



(57) 



Equations (56) can be integrated to give the ray trajec- 
tory 



dx 



±^/ S m 2 (m/iz) 2 - 1, 



where I and l z are the constants of integration in ( 56 1 , and 
the sign of the derivative changes when the square root 
goes to zero. Noting the symmetry of the medium, and 
setting $ = 7r/2, I = l z , it is clear that rays propagate 
along the great circles and polarization has a constant 
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orientation along each ray (figure [2] a) . In Riemannian 
transformation optics, we can only manipulate the prop- 
agation of polarization along a ray through changing the 
geodesies. However, the additional geometrical property 
of torsion, outlined in sections |IV| and |VB| allows us to 
'twist' the co-ordinate lines on the sphere so that the 
polarization is changed, while leaving the geodesies un- 
altered. 

For instance, if we wished to generate TE polarization 
at one point on the sphere (e.g. = tt/2, $ = tt/2, x = 
0), and have it arrive at the antipode (9 = 7r/2, <E> = 
ir/2,x — 7r) with TM polarization, then we could add 
the following torsion into the connection (this picture is 
accurate so long as a » A), 



T, 



'jf 



2 



(58) 



The behaviour of the polarization on the surface of the 
2-sphere with the addition of ( 58 1 is shown in figure [2j 



inhomogeneous media with isotropic chirality. If the 
definition of the field tensor is kept in the same gauge 
invariant form as in a Riemannian geometry, then we 
have shown that a background with a particular kind of 
torsion can reproduce the correct constitutive relations 
for such media. Furthermore, if the chirality produces 
only a small amount of optical activity over a single 
optical cycle, and varies by only a small amount over a 
wavelength, then we have also shown that geometrical 
optics behaves exactly as if on a background with a 
non-zero torsion. This formalism allows for a greater 
degree of control over polarization within the theory of 
transformation optics, and we have given an example 
of how torsion can be implemented to the change of 
polarization without changing the geodesies. 
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